www.gusucode.com > 基于matlab编程蒙特卡罗仿真设计源码程序 > 基于matlab编程蒙特卡罗仿真设计源码程序/第五讲 蒙特卡洛方法的应用/litiR2.m

    function Rguji=litiR2(t,thetaa1,thetaa2,thetab1,thetab2,mm)
%t 是要求统存寿
%thetaa1 是元件A1的数学期要
%thetaa2 是元件A2的数学期要
%thetab1 是元件B1的数学期要
%thetab2 是元件B2的数学期要
%mm 是随机实验次数
frq=0;
randnuma1 = exprnd(thetaa1,1,mm);
randnuma2 = exprnd(thetaa2,1,mm);
randnumb1 = exprnd(thetab1,1,mm);
randnumb2 = exprnd(thetab2,1,mm);
for ii=1:mm
    if (randnuma1(1,ii)>t)|(randnumb1(1,ii)>t)
        pass1=1;
    else
        pass1=0;
    end 
    if (randnuma2(1,ii)>t)|(randnumb2(1,ii)>t)
        pass2=1;
    else
        pass2=0;
    end 
    if (pass1*pass2)==1
        frq=frq+1;
    end    
end    
Rguji=frq/mm